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A recently proposed two orbital model for the new Fe-based superconductors is studied using the 
Lanczos method on small clusters as well as pairing mean-field approximations. Our main goals 
are (i) to provide a comprehensive analysis of this model using numerical techniques with focus on 
the magnetic state at half-filling and the quantum numbers of the state with two more electrons 
than half- filling and (ii) to investigate the nodal structure of the mean-field superconducting state 
and compare the results with angle-resolved photoemission data. In particular, we provide evidence 
that the dominant magnetic state at half-filling contains spin "stripes" , as observed experimentally 
using neutron scattering techniques. Competing spin states are also investigated. The symmetry 
properties of the state with two more electrons added to half filling are also studied; depending 
on parameters, either a spin singlet or spin triplet state is obtained. Since experiments suggest 
spin singlet pairs, our focus is on this state. Under rotations, the spin-singlet state transforms as 
the _B2g representation of the D^i, group. We also show that the s± pairing operator transforms 
according to the ^ig representation of Diy, and becomes dominant only in an unphysical regime of 
the model where the undoped state is an insulator. We obtain qualitatively very similar results both 
with hopping amplitudes derived from a Slater-Koster approximation and with hoppings selected 
to fit band-structure calculations, the main difference between the two being the size of the Fermi 
surface pockets. For robust values of the effective electronic attraction producing the Cooper pairs, 
assumption compatible with recent angle-resolved photoemission (ARPES) results that suggest a 
small Cooper-pair size, the nodes of the two-orbital model are found to be located only at the 
electron pockets. Note that recent ARPES efforts have searched for nodes at the hole pockets or 
only in a few directions at the electron pockets. Thus, our results for the nodal distribution will 
help to guide future ARPES experiments in their search for the existence of nodes in the new Fe- 
based superconductors. More in general, the investigations reported here aim to establish several of 
the properties of the two orbital model. Only a detailed comparison with experiments will clarify 
whether this simple model is or not a good approximation to describe the Fe pnictides. 



PACS numbers: 74.20.Mn,74.20.-z,71.27.+a 

I. INTRODUCTION 

1. Current status of experimental and theoretical 
investigations 

The discovery of a new family of superconducting ma- 
terials with Fe-As layers in their structureii^i^i^i^i^iii^ has 
triggered a large effort in the condensed matter com- 
munity. LaOi _2;F3;FeAs is a much studied example of 
this family of compounds. The ~ 55 K record critical 
temperature^ in SmOj _3:F2;FeAs is second only to those 
observed in the Cu-oxide family of high temperature su- 
perconductors. In addition, there are several aspects of 
the physics of the new Fe-based superconductors that 
suggest the possibility of an exotic pairing mechanism at 
work: 

(1) Evidence is accumulating that phonons may not be 
sufficient to understand the superconductivity of these 
compounds 1^ 1^°'^^ Moreover, the importance of correla- 
tions between the electrons has been remarked in sev- 
eral investigationsffi^^i^^iii^i^^^^ In fact, it 
has been claimed that these oxypnictide superconduc- 



tors may bridge the gap between MgB2 and the Cu- 
oxide superconductors i^^i^ In addition, a pseudogap was 
detected, similarly as in the cupratesj ^^'^^i^^'^^ Coexis- 
tence or proximity of magnetism and superconductivity 
has also been report edj ^^i'^°'^^i'^^ Although the parent un- 
doped compound is not a Mott insulator, these results 
suggest that the influence of electron-electron repulsions 
cannot be neglected. Perhaps the intermediate range of 
where U is the typical Hubbard repulsion scale 
and t the typical hopping amplitude in a tight-binding 
description, is the most representative of the new super- 
conductors. U cannot be too large, otherwise the system 
would develop a gap and the undoped compound would 
be insulating, contrary to the experimentally observed 
properties of the undoped limit that suggest bad metal- 
lic behavior. But poor-metal characteristics imply that U 
cannot be too small either, otherwise the undoped system 
would be a good metal. In addition, the mere presence 
of a spin-density-wave magnetic state shows that corre- 
lations must be important. 

(2) Several experimental investigations suggest 
the presence of nodes in the superconducting 
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gap^^^^^^^^d2^ This is reminiscent of the 
nodes that appear in the d-wave superconducting state 
of the high-Tc cuprates. However, other investigations 
indicate nodeless superconductivity] As 
a consequence, this issue is still controversial. 

(3) The undoped parent compound has long-range spin 
order in the ground statci^ This magnetic state corre- 
sponds to spin "stripes" having the Fe spins along one of 
the Fe-Fe crystal axes pointing all in the same direction, 
and being antiferromagnetically coupled in the perpen- 
dicular direction. According to neutron scattering exper- 
iments, in LaOFeAs the transition to this magnetic state 
occurs at 134 K, and the magnetic moment is 0.36 ^b, 
which is smaller than anticipatedi^ For NdOFeAs^^ the 
critical temperature is 141 K and the magnetic moment 
is even smaller 0.25 ^b- However, recently by means of 
resistivity, specific heat, and magnetic susceptibility mea- 
surements, the antiferromagnetic critical temperature of 
SrFc2As2 was reported to be as high as 205 K, with a 
more robust Fe magnetic moment of value 1.7 ^b*^ Also, 
CaFe2As2 was investigated using neutron diffraction, and 
a critical temperature 173 K with a moment 0.8 fXB was 
reportedi^ Thus, although originally it was believed that 
the undoped material had a very weak magnetic state, 
the most recent results suggest that the spin striped or- 
der may be more robust. 

On the theory front, several band-structure calcula- 
tions have shown that the Fermi surface of these and 
related compounds is made out of two small hole pock- 
ets centered at the F point, and small electron pockets 
at the X and Y points, in the notation corresponding to 
a square lattice of Fe atomsi ^^i^'''^^i^^'^° These calcula- 
tions have also shown that the 3d levels of Fe play the 
dominant role in establishing the properties of these ma- 
terials near the Fermi level. To address theoretically the 
physics of these compounds, particularly the supercon- 
ducting state, model Hamiltonians arc needed and sev- 
eral proposals for the dominant pairing tendencies have 
been made^'^'S'^*^'^^^^^^^ In particular, a 
two orbital model based on the dxz and dyz orbitals was 
recently presented.— Several other investigations have 
addressed this model for the new superconductors, using 
a variety of approximations ^^'Ze^'Z^S'^^^^^^ 
Classifications of the possible superconducting or- 
der parameters for the two-orbital model have been 
madeMi^s^ 

As already mentioned, a variety of experimental results 
suggest that the Cooper pairs are spin singletsj ^^'^°i^° 
Thus, it is important to find the range of parameters lead- 
ing to spin singlets in model Hamiltonians, since several 
calculations produce either singlet or triplet supercon- 
ductivity depending on the couplings and bandwidths 
used. For this experimentally-based reason, our focus 
here will be mainly on singlet superconductivity. 



2. Why the two orbital model? 

In this manuscript, a detailed study of the two orbital 
model for the oxypnictide superconductors is carried out 
using Lanczos and pairing mean- field techniques. This 
effort provides a comprehensive view of the model, con- 
siderably expanding our recent research on the subject^^ 
by varying the several couplings of the model and study- 
ing the main tendencies. When two electrons are added 
to the half-filled ground state, a spin-singlet state that 
transforms in a non-trivial manner under rotations is 
shown to dominate in the regime of couplings that is 
argued to be the most relevant to describe the new super- 
conductors. In addition, the nodal structure of the super- 
conducting state obtained using these spin-singlet pairs is 
here studied for this model using the pairing mean-field 
approximation. Our results are compared with recent 
ARPES experiments, and suggestions to further refine 
the search for nodes in those experiments are discussed. 

Currently there is no consensus on what is the min- 
imal model capable of capturing the essential physics 
of the oxypnictides. Band structure calculations in 
the local-density approximation (LDA) indicate that the 
bands that form the observed electron and hole pock- 
ets are strongly hybridized but they have mostly FeSd 
character Several authors argue that the hybridiza- 
tion of the Fe-3(i is so strong that all 5 d orbitals have 
to be considered to construct a minimal model. For in- 
stance, a five-orbital model has been proposedi^ The 
tight-binding term respects the FeAs lattice symmetries 
and the hopping parameters have been obtained from 
fittings against the LDA calculations. The parameters 
used reproduce the Fermi surface (FS) for the electron 
doped system (i.e. electronic density n = 6.1) but an 
extra hole pocket around M (in the notation of the ex- 
tended Brillouin zone) appears for the undoped case and 
upon hole doping. For this reason, the model may not be 
suitable to study the magnetic properties of the undoped 
system. In addition, the number of degrees of freedom 
in fivc-orbital models makes its study very difficult using 
numerical techniques. However, LDA calculations have 
shown that, although heavily hybridized, the main char- 
acter of the bands that determine the FS is dxz and dyz, 
with a small contribution of d^y at the most elongated 
portions of the electron pockets i^SiSi This fact has been 
the main justification for the proposal of twoSiZ^ and 
three^ orbital models. The two-orbital model can have 
its hopping parameters fitted such that the shape of the 
FS, both in the undoped and electron and hole doped 
cases, arc well reproduced in the reduced or folded BZ. 
However, it has been argued by some authors^ that the 
two hole-pockets around F have to arise from the dxz and 
dyz orbitals that are degenerate with each other at F, as 
obtained in LDA. In the two-orbital model, one of the 
hole pockets forms around M in the extended BZ which 
gets mapped onto F upon folding. The d^z and dyz or- 
bitals that form the M-point pocket are degenerate at M 
and, upon the folding, give rise to higher energy bands 



3 



at the r point. For this reason one of the hole pockets 
in the two orbital model may not have the correct linear 
combination of orbitals, potentially leading to incorrect 
conclusions. In addition, it is also argued that the contri- 
bution of the dx_y orbital to the electron pockets may play 
an important role that should not be ignored which moti- 
vated the proposal of the three-orbital modeli^ However, 
the three orbital model cannot eliminate a spurious hole 
pocket around M. Thus, a fourth orbital needs to be 
added to accomplish this task and, again, the number of 
degrees of freedom makes this model too complex to be 
studied numerically. 

Then, the justification for continuing studying a min- 
imal model with just two orbitals, as carried out in the 
present manuscript, is the following: (i) The correct 
shape of the FS is reproduced in the reduced Brillouin 
zone, both in the doped and undoped cases, (ii) The 
main character of all the bands that determine the FS 
is dxz and dyz, except for a small portion of the electron 
pockets that has dxy character. Then, it it worthwhile 
to understand the role, if any, that this orbital plays in 
the magnetic and superconducting states, (in) The two- 
orbital model is the only one that can be studied exactly 
with numerical techniques using the minimal size cluster 
needed for a spin striped state.— Thus, we believe that 
it is very important to establish which properties of the 
oxypnictides are properly captured by this model, and 
which ones not. The role that the correct shape of the 
FS plays can be investigated as well, and also the pair- 
ing symmetry and nodal structure involving only the d^z 
and dyz orbitals. It is interesting to notice that although 
two superconducting gaps may appear in a two orbital 
model,— symmetry forces the magnitude of the gaps to 
be the same in this case<^ 



3. Organization 

The organization of the paper is as follows. In Sec- 
tions II and III, the two-orbital model is derived. The 
emphasis is on the Slater-Koster (SK) procedure to eval- 
uate the hopping amplitudes, but the model derived by 
this method is more general: it coincides with the two- 
orbital Hamiltonian proposed earlier,— and the values 
of the hoppings can be obtained also by fitting band- 
structure calculationsi^^ Both sets of hopping parameters 
will be used in the following sections. The qualitative as- 
pects of the magnetic and pairing states are shown to be 
the same for both sets of hopping amplitudes. In Section 
IV, results for the ground states of the undoped model 
(half-filled) and the case of two more electrons than half- 
filling will be discussed using the Lanczos technique. The 
emphasis is on the dominant magnetic states and on the 
pairing tendencies, which are either in the spin singlet 
or triplet channels depending on couplings. Moreover, 
the spin singlet case is shown to correspond to the B2g 
representation of the D^ii lattice symmetry group of the 
model. Section V contains a pairing mean-field analysis 



of the nodal structure of the model. The two orbitals 
nature of the problem causes the number and location 
of the nodes to be a more complex topic than for just 
one orbital. A qualitative comparison with experiments 
is included here. Section VI contains our main conclu- 
sions. The possible source of the i?2g pairing and the s± 
pairing operator are discussed in the appendices. 



II. MODEL DISCUSSION AND DERIVATION 
OF HOPPING AMPLITUDES 

To study numerically the properties of LaO j _:EF:EFeAs 
and related compounds, it is necessary to construct a 
simple model, one that contains a minimum amount of 
degrees of freedom but still preserves the main physics 
of the problem. Since all the materials in the family 
have in common the Fe — X planes (X=As, P, ...), as a 
first approximation we will just focus on those planes, 
similarly as it occurs in theoretical studies of the Cu- 
O planes in the cuprate superconductors. In addition, 
band structure calculation o^^i^^'^^i^^'^° have shown the 
relevance of the Fe 3d levels, and that mainly two bands 
determine the Fermi surface (see Introduction). Based 
on these considerations, here we will include only the dxz 
and dyz Fe orbitals in our discussion. To estimate the 
hopping amplitudes for the tunneling from one Fe to an- 
other and, thus, define a tight-binding model, we will cal- 
culate their hybridization with the three p orbitals of As 
following the Slater-Koster formalism.— From the Fe-As 
hopping integrals, we will calculate the effective Fe-Fe 
tight-binding hopping parameters following a standard 
perturbative approach. Thus, the hopping parameters in 
this model will be functions of the overlap integrals be- 
tween the orbitals and the distance between the atoms. 
While this procedure is not as accurate as band-structure 
calculations, it provides a simple to understand approach 
that has "ab-initio" characteristics, can be easily repro- 
duced since the calculations are analytical, and they also 
illustrate how the geometry of the problem affects the 
hoppings. 

However, before proceeding, we remark that another 
avenue to obtain the hopping amplitudes is via fittings 
of the band-structure calculations.— In our description of 
results below, data for both the SK hoppings and those 
that fit band structures will be presented. An important 
result is that both sets of hoppings lead to similar qual- 
itative results, both in the undoped case, regarding the 
magnetic state, as for two electrons added, regarding the 
pairing tendencies. 

The unit cell in the FeAs planes contains two Fe atoms, 
since the As atoms are above and below the plane defined 
by the Fe atoms in alternating plaquettes (Fig. [ija)). 
However, after the calculation previously described only 
the Fe atoms will be considered in a simple two-orbital 
Hamiltonian. Since these Fe atoms form a planar square 
lattice, it is natural to orient the lattice as in Fig.[T]^b). 

To guide the discussion, consider a cluster with 4 Fe 
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TABLE I: Coordinates of 
the atoms in Fig. [21 a). 
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TABLE II: Director 
cosines of the Fe atoms 
with respect to Aso in 
Fig. Ha). 



FIG. 1: (Color online) (a) Schematic representation of the Fe- 
As plane. Blue circles are the Fe atoms. The red filled circle is 
an As atom at a distance c below the plane, while the red open 
circles are As atoms at a distance c above the plane, (b) Unit 
cell for the effective Fe-only square lattice. The Fe-Fe lattice 
has been rotated by 45°. (c) Schematic first Brillouin Zone 
(FBZ) for the Fe-As plane. The point X is at (27r/d, 0), with 
d = -spil. (d) FBZ for the Fe-As lattice after a 45° rotation, 
(e) FBZ for the rotated Fe-Fe shown in (b). X' = (27r/Z,0) 
and it is equivalent to the M point for the Fe-As plane in 
(c). The electron and hole Fermi surfaces obtained by band- 
structure calculations are schematically indicated. Panels c-e 
will be useful for the discussion related to the nodal structure 
of the superconducting state in Section V. 



and 5 As atoms (Fig.[21^a)). The coordinates of the atoms 
are needed to calculate hopping amplitudes, and they are 
provided in Table I, where fc, and c are obtained from 
the materials structure. The nearest-neighbor (NN) Fe- 
Fe distance is / = 2.854 thus k = 1/2 ^ 1.427 A. 
The distance between Fe and As is s = 2.327 AjS see 
Fig. H^b). The next-nearest-neighbor (NNN) Fe-Fe dis- 
tance along the square diagonal is d = \/2l = 4.037 A(see 
Vig.^c)) and r = d/2 = 2.018 A. According to Fig.[2td), 

1.158 A. The director 



- F/2 



cosines I, m, and n for each of the Fe atoms^S^ with re- 
spect to the As located at (0,0,-c), are given in Table IT 



A. Overlap integrals between the Fe dxz and dyz 
orbitals and the As and Py orbitals 



(a) 



(b) 



As 



O As, 



Fe 



Fe, 



Fe. 



As„ As, 

^ ^ 



Fe, 



OAs4 



Fe. 



Fe 



(d) 



d.-' 



Fe. 



Fe, 



As 



Fe 



r=d/2 



FIG. 2: (Color online) (a) The Fe-As cluster used in our cal- 
culations of the hoppings. Green circles are the Fe atoms. 
The red circle at the center is an As atom at a distance c be- 
low the plane, while the shaded red circles are the As atoms 
that are a distance c above the plane, (b) Distances s and I 
for NN Fe-Fe atoms, (c) The distance d along the diagonal of 
the Fe-Fe plaquettes. (d) The distance c for As atoms. 



According to the SK analysis, for the orbitals consid- 
ered here we obtain the following results for the center 
integrals: 
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V3lmn{pda) — 2lmn{pd7r) , 
\/3lmn{pda) — 2lmn{pd'!r) , 
V3l^n{pda) + n{l - 2f){pd'K), 

(1 - 2m'^)(pdTT). 
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(3) 
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The corresponding hopping amplitudes are 



\tx.-. 



\ty.xz I 



^/2,^{pda) - 2^{pdTr), (5) 



\tx,xz\ = \ty,yz\ - b = Vi^ipda) + ^(1 - 2— )(pd7r).(6) 

The signs and values of these hoppings for the cluster 
that we are considering are in Fig. 3. The values of I, to, 
and n shown in Table I are for Asg , while some signs will 
be different for Asi, As2, Asa, and AS4. 

Using the values of fc, s, and c given above, we obtain: 



a = 0.324(pdcr) - 0.374(pd7r), 
b = 0.324(pdcr) + 0.123(pd7r). 



(7) 
(8) 
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FIG. 3: (Color online) (a) Hoppings between dyz (dxz) orbitals 
in Fe and Px (Py) orbitals in As for the cluster considered in 
Fig-Oa). (b) Hoppings between dyz (dxz) orbitals in Fe and 
Py (px) orbitals in As for the cluster considered in Fig[2ja). 
Continuous (dashed) lines indicate positive (negative) values. 



the inter-orbital hopping. The fact that the plaquette 
diagonals are equivalent by symmetry implies that the 
absolute values of the hoppings must be the same along 
these diagonals, but the signs can be different as shown 
here. More exphcitly, we obtain: t^^^^ = ab/A, and 
-ab/A. 



,xz-yz _ 



B. Overlap between dxz and dyz with pz 

The consideration of the orbitals adds two new cen- 
ter integrals to the present analysis: 

Ez,xz = Vin\pd(j) + l{l-2r?){pdTi), (10) 
Ez.yz = VSn^m(pda) + m{l — 2n^){pdTT), (11) 

which means that a new hopping must be considered 



Now let us compute the hopping amplitudes for a 
square lattice made up only of Fe atoms. For the NN 
effective Fe-Fc hopping i„„ we will consider the pair of 
atoms Fei and Fe2. For the hopping between the dxz or- 
bitals, there are two possible paths using the As Py as a 
bridge. Their contribution is given by (1) dxzFei-pyAso- 
dxzFc2 and (2) (ixzFei-pyAsi-dxzFG2. From Fig. [3]Ja), we 
observe that these paths contribute with — each. Re- 
garding the use of the Px of As as a bridge, in this case 
there are also two paths: (3) c?xzFei-pxAso-dxzFe2 and (4) 
dxzFei-pxAsi-dxzFe2. From Fig. O^b), these paths con- 
tribute with b^ each. Reasoning in an analogous man- 
ner, four similar paths are found for the NN hopping 
between orbitals dyz: (1) dyzFei-pxAso-dyzFe2 and (2) 
dyzFei-pxAsi-dyzFe2, that from Fig.[3]Ja) they give a con- 
tribution —a^ each, and (3) dyzFei-pyAso-dyzFe2 and (4) 
dyzFei-pyAsi-c?yzFe2, that from Fig.[21Ib) they give a con- 
tribution b^ each. Combining all these results, and to 
second order in perturbation theory^ the Fe-Fe nearest- 
neighbor hopping amplitude is given by: 



(-2a2 + 262)/A = 2(62_a2)/A, (9) 



where A is the difference between the on-site energies 
of the d and p orbitals. Notice that by mere geometri- 
cal reasons, it is not possible to have a nearest-neighbor 
hopping from dyz to dxz. 

For the hopping td along the Fe lattice plaquette diag- 
onal, namely the NNN Fe-Fe hopping, let us consider the 
hopping from Fci to Fcs and from Fe2 to FC4. It can be 
easily shown that dxzFei-pxAsg-dxzFcs contributes by an 
amount 6^ to t^^, while dxzFci-pyAso-dxzFea contributes 

to . The same result is obtained if the hopping from 
Fe2 to Fe4 is considered. Combining these numbers, then 
we obtain = tf = (a^ + b^)/A. 

Along the plaquette diagonal we can also obtain inter- 
orbital hopping. From Fci to Fes the contribution is -a6, 
while from Fe2 to Fe4 it is ab. Thus, the hopping along 
the x + y and x — y directions are different by a sign from 



(12) 



\tz.xz\ = \tz,yz\ =g = V3-^{pda) + -(1 - 2^)(pd7r). 
Using the values for fc, s, and c calculated before, 

g = 0.263(pd(T) + 0.31(pd7r). (13) 
The signs are indicated in Figs.|4t^a) and (b). 
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FIG. 4: (Color online) (a) Hoppings between the dyz orbitals 
in Fe and the pz orbitals in As, for the cluster considered in 
Fig. [2ja) . (b) Hoppings between the dxz orbitals in Fe and 
the Pz orbitals in As, for the cluster considered in Fig. [2Ia)- 
Continuous (dashed) lines indicate positive (negative) values. 



Thus, we obtain an additional contribution to the NN 
hopping i„„ so that = Ig^ j Al {-2g'^ / A') along the 
y {x) axis. Reciprocally, = ^g"^ / A' {~2g^/A') along 
the X (y) axis. Along the diagonal, td ~ —g^/A' for 
both orbitals is obtained. Note also that Pz generates an 
inter-orbital diagonal hopping given by —g^jA! (t/^/A') 
along the x ^ y (x — y) directions. A' is the difference 
between the on-site energies of the d and Pz orbitals. 
From Ref. [H, the gaps are A = 1.25 eV and A' = 5 eV, 
but other values for these gaps are also considered below. 
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C. Direct Fe-Fe hopping 

Since the distance between Fe atoms is I = 2.854 A, 
comparable to the Fe-As distance, the contributions to 
the electron hoppings coming from the direct overlap be- 
tween the d orbitals of the Fe atoms should also be con- 



sidered. Following SK,— E. 



APn^){ddT:) + (m^ + Pn^)ddS, Ey.^y, = 3m^n^{dda) + 
{m^ ~ Am'^n^){ddT:) + {l^ + m?n^)dd5, and Exz,yz = 
3lmn^(dda) + Im (1 - 47i2) [(ddTr) - {dd5)] . Notice that all 
the Fe atoms have n = 0, and I — ±1, m = (Z = 0, 
m = ±1) if they are neighbors along the x (y) direction. 
Thus, the inter-orbital hopping vanishes, and we obtain 
txz,xz = —ddiT {tyz,yz ~ —ddn) along the direction x 
(y), and ddS along y (x). These same expressions can be 
used to obtain the diagonal Fe-Fe hopping parameters. 
We find that q, = —{ddir' + dd5')/2, where a = xz or 
yz, while t'^z,yz = ^iddir' — ddS')/2 with the minus (plus) 
sign for the x + y (x — y) direction and the prime indicates 
second nearest- neighbors overlap integrals. 



III. EFFECTIVE TWO-ORBITAL 
TIGHT-BINDING MODEL 



A. Hopping Term 

Considering the results of the previous section, the 
kinetic-energy term of the effective tight-binding Hamil- 
tonian involving the d^z and dyz orbitals, defined on the 
square lattice formed only by the Fe atoms, is given by: 

HtB = -hJ2(4,x,adi+y,x,a + dly^^di+i^y^a + h.c.) 



~t3 ^ {dl,^^^di+fi+C'.x,a + dly^^di+fX+p^y^^ + h.c. 



'ti'^i4,x,cr'^i+x+y,y,<T + dly^^di+x+y,x,a + h.C.) 



ti ^~^('^i.x.a-'^'+^-y.y.g '4,y,a^i+x-y.xM + h.C. 



(14) 



In this Hamiltonian, the operator dl ^ ^ creates an elec- 
tron with spin z-axis projection cr, orbital a, and on the 
site i of a square lattice. The chemical potential is given 



by fi and nf are number operators. The index fi = x 
or y is a unit vector linking nearest-neighbor sites. The 
hoppings, within the SK approach, are given by: 



-2[(62 -a2)/A + .gVA'] -dd<5, 

-2[(62 - a2)/A - .gVA'] - ddii, 

-[(a^ + b'^A - g^A')] - {ddir' + ddS')/2, 

-{ab/A - g^/A') - {ddir' - dd6')/2. (15) 



The explicit expressions for these hopping amplitudes in 
terms of the overlap integrals using the parameters for 
FeAs can be easily found and they will not be provided 
here. The two orbital model proposed by Raghu et alM 
has the same form as the one presented above but the 
hoppings are obtained by fitting band structures<^ 

It is interesting to notice that if only the direct over- 
lap between the d orbitals is considered, i.e. ignoring the 
indirect hopping through the p As orbitals, the form of 
Eq. does not change. Thus, the form of 7?tb arises 
from the symmetry properties of the dxz and dyz orbitals 
rather than from the location of the As ions. However, 
the indirect Fe-Fe hopping through the As atoms plays 
a key role in providing the relatively large value of the 
diagonal hopping ^3 vs. the NN hoppings which, as dis- 
cussed below, stabilizes the magnetic stripe order. For 
example, if we only consider the direct hopping then 
^3/^2 ~ ddir' /2ddTT, where dd6' w was assumedi^ Since 
ddn' <C ddir, then jtaj ^ \t2\. However, if wc consider 
the indirect hopping then jtsj > |t2|/2. 




FIG. 5; (Color online) (a) Energy vs. momentum for the 
non-interacting tight-binding Hamiltonian in Eq. (|14|l us- 
ing ti = 0.058 eV, t2 = 0.22 eV, ts = -0.21 eV, and 
t4 = —0.08 eV. These hopping amplitudes are obtained from 
the Slater-Koster formulas using (pda) = 1 eV and (pdn) = 
—0.2 eV, supplemented by A = A' = 1 eV for simplicity, and 
H — —0.03 eV, which corresponds to half-filled orbitals. Re- 
sults are plotted along the path (0, 0) — (vr, 0) — (tt, tt) — (0, 0). 
(b) Fermi surface for the half-filled system. 



To analyze the influence of the several parameters, let 
us consider two special cases. Setting (pda) = 1.0 eV, 



7 



pdTT = -0.2 eV, A = 1.0 eV, and A'=l eV in Eqs. (IS- 
IS), and neglecting the direct Fe-Fe coupling i.e. using 
ddTT = ddS = 0, we obtain: h = 0.058 eV, t2 = 0.22 eV, 
h = -0.21 eV, and i4 = -0.08 eV. With these values, the 
band structure, shown in Fig.[5l is qualitatively similar to 
the band-structure calculations, although the pockets are 
larger in size. Another example can be obtained by using 
the calculated values of the energy gaps, which arc A = 
1.25 eV and A' = 5 eV,^ In Fig. [H the band structure 
is shown for {pda) = 1 eV, pdn = -0.2 eV, A = 1.25 eV, 
A' = 5 eV, dd-K = 0.2 eV, and ddS = -0.02 eV, i.e. 
including the direct Fe-Fe hopping. Now the hole pocket 
at r is larger than in the previous case, but the overall 
shape remains similar. 
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FIG. 6: (Color online) (a) Energy vs. momentum for 
the non-interacting tight-binding Hamiltonian Eq. (|14|l us- 
ing ti = -0.1051 eV, t2 = 0.1472 eV, ts = -0.1909 eV, and 
ti — —0.0874 eV, obtained using the parameters (pda) — 
1 eV, (pdiv) = -0.2 eV, ddn = 0.2 eV, ddS = -0.02 eV, 
A = 1.25 eV, and A' = 5 eV. The chemical potential is 
^ — 0.081 eV, which corresponds to half-filled orbitals. Re- 
sults are plotted along the path (0, 0)-(7r, 0)- (tt, 7r)-(0, 0). (b) 
Fermi surface for the half-filled system. 



Notice that the overlap integrals can also be estimated 
using tabulated values and the distances between the 
atoms 1^ The band structure and Fermi surface obtained 
using these values are shown in Fig. [T] 

To complete the analysis, let us discuss now the results 
obtained using the set of hoppings that fit band-structure 
calculations The dispersion and Fermi surface arc in 
Fig. [8l By construction, the agreement with the band- 
structure Fermi surface is better than in the other cases, 
the main difference being the size of the pockets. Never- 
theless, it appears that in a broad range of hoppings and 
couplings, the qualitative topology of the Fermi surfaces 
remains the same, and this is probably the reason why 
the main magnetic and pairing properties are also similar 
among the many sets, as shown explicitly below. 




FIG. 7; (Color online) (a) Energy vs. momentum for the 
non-interacting tight-binding Hamiltonian in Eq. (|14|l using 
ti = -0.129 eV, t2 = 0.05 eV, ts = -0.137 eV, and t4 = 
-0.019 eV obtained for (pda) = -0.41 eV, (pdvr) = 0.19 eV, 
ddn = 0.18 eV, ddS = 0, A = 1.25 eV, and A' = 5 eV. The 
chemical potential is ^ = 0.053 eV, which corresponds to half- 
filled orbitals. Results are plotted along the path (0, 0)-(7r, 0)- 
(tt, 7r)-(0, 0). (b) Fermi surface for the half-filled system. 




FIG. 8: (Color online) (a) Energy vs. momentum for 
the non-interacting tight-binding Hamiltonian Eq. (|14|l using 
the hopping amplitudes obtained from fits of band-structure 
calculations'^ ti = —1.0, t2 = 1.3, and = ti = —0.85 (all 
in eV units). Results are plotted along the path (0, 0)-(7r, 0)- 
(tt, 7r)-(0, 0). (b) Fermi surface for the half-filled system. 



B. Interactions 

In this section, the Coulombic interaction terms are 
added to the tight-binding Hamiltonian Eq. (fT4)) to form 
the full two-orbital model. These Coulombic terms arei^ 

Hint = "-i,a,T"-i,aa + ^ J/"^) X! "i.^"'',!' 



2J Su^-Si^y+J Yidl^^'fdl^^^di,y,idi,y,'\+h.c.). (16) 
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where a = x,y denotes the orbital, Si^Q {'^1,0) is the spin 
(electronic density) in orbital a at site i, and we have used 
the relation U' = U — 2 J, from rotational invariance.— 



IV. EXACT DIAGONALIZATION RESULTS 



Method 



C. Pairing 

Diagonalizing exactly the full two-orbital Hamiltonian 
on a X V8 cluster with periodic boundary conditions, 
it was observed in previous investigations that in regions 
of parameter space the ground state with two extra elec- 
trons above half-filling is a spin tripleti^ In this case, the 
relevant pairing operator is given by: 

At(i). = Y.(<^A+,,y.^ - iyA+,,^J^ (17) 



where a =t or | denotes the spin projection 1 or -1, 
respectively, while the projection operator is: 



4,y,A+,,^,i 4.a^H-M,..T)' (18) 

or, in momentum space, 

At(k), = (cosfcx+cosfc,y)(4^^^dlk,^,^-4^^^dt_^,^ 

(19) 



In this section, the Lanczos or Exact Diagonalization 
(ED) method will be used to obtain the ground state of 
the two-orbital model, both at half filling and also for a 
system with two electrons more than half filling. Due to 
the exponential growth of the Hilbert space with increas- 
ing cluster sizes, here our effort must be restricted to a 
tilted -x/Sxa/S cluster <^ Using translational invariance, 
the Hilbert space can be reduced to 21, 081, 060 states at 
half filling and 16, 359, 200 for two electrons away from 
half filling. Taking into account the additional symme- 
tries of spin inversion as well as rotations, the dimension 
of the Hilbert space becomes 2, 600, 000. The employed 
Lanczos scheme is standard and requires up to 11 GB of 
memory when only translational invariance is used. Note 
that the two-orbital 8-sites cluster has a similar Hilbert- 
space size as a 16 sites one-band Hubbard lattice, and 
they are similarly computationally demanding. The fo- 
cus of our effort is on ground states for a fixed set of quan- 
tum numbers corresponding to the symmetries that were 
implemented. We use both the hoppings from the SK 
approach and also the hoppings that fit band-structure 
calculations, and find qualitatively consistent results for 
both sets. 



A'''(k)o = (cos/cj; + cos ky)(dl^^-^d\y^^ + 

<..i^-k,,,T - 4,,,T^-k,..i - dly,idU.^,l)- (20) 

This operator is invariant under the irreducible rep- 
resentation of the group D^h, it is odd under orbital ex- 
change, and it is a spin triplet. 

However, in our previous effort we have also identified 
regions of parameter space where the state with two extra 
electrons is a spin singlet, which appears to be compatible 
with the results of experiments that favor singlet states 
over triplets.— i^i^^ The dominant pairing operator for 
the singlet is given by: 



d 



that in momentum space becomes 



A'''(k) = ^(cosfcj; +cosky)dl^^^dl^_ 



(21) 



(22) 



This operator transforms as the B2g irreducible repre- 
sentation of the point group, it is even under orbital 
exchange, and it is a spin singlet. As explained in the 
introduction, the experimental results favoring spin sin- 
glet pairing lead us to focus our effort on this spin-singlet 
operator in the following sections. 



B. Results using Slater-Koster derived hoppings 

1. Fermi surfaces, spin order, and spin of the pairs 

In the SK approach, the parameter pdcr is here kept 
fixed equal to 1, providing the scale, and the free param- 
eter pdn is varied. To constrain the values of pdn, let us 
return to the tight-binding Hamiltonian. Figure [H] shows 
how the Fermi surface evolves by changing pd?:. These 
figures are in the unfolded Brillouin zone, i.e., for one Fe 
atom per unit cell. For a negative ratio pdir/pda, hole 
pockets around momenta (0, 0) and (tt, tt) and electron 
pockets around (0, 7r)/(7r, 0) are found. However, for a 
vanishing ratio pdi: = 0, additional electron pockets ap- 
pear at (7r/2,7r/2), while the pockets around (0, tt) and 
(tt, 0) disappear fast by further increasing pdn to posi- 
tive values. As discussed before,— the robust NNN hop- 
ping t3 at negative pdir/pda induces tendencies toward 
a (0, tt) / (tt, 0) magnetic ordering at half filling, as shown 
by the spin structure factor in Fig. 



10(a) 



This is in good 
agreement with neutron scattering experiments. Thus, it 
is clear that the realistic regime corresponds to negative 
pdir, and an opposite sign of pdi: and pda is also what 
would be expected from the tabulated values 1^ 

As it can be observed in Fig. 
U enhances the spin "striped" 
dominant even at U= 



10(a) 



the onsite repulsion 
ordering, which is already 
=0 (although in this nonintcracting 



case a power-law decay in the spin correlations is ex- 



9 



(a) pdn/pda = —0.5 
1r 



> 




(b) pdir/pda = -0.2 
1r 

0.5 



-0.5 
-1 




-1 





k /7t 



0.5 



--^1 

(e) 
1 


1 

k /7i 

X 

pdir/pda = 0.1 


0.5 















-0.5 












k /7i 

X 


1 



(d) pdn = 
5 • • 

0.5 • • 

h r 







k /jt 

X 

(f) pdn/pda = 0.5 
1 



0.5 



-0.5 



1 




1 




f 










J 




1 




r 
1 



k /7t 



FIG. 9: (Color online) Non-interacting {U = J = 0) Fermi 
surface in the unfolded Brillouin zone, at the pdn/pdas indi- 
cated. The realistic regime is pdn/pda < since hole pockets 
around (0,0) and (7r,7r), and electron pockets around (0, tt) 
and (0,7r), are observed. 
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FIG. 10: (Color online) Spin structure factor S{k) at half 
filling obtained with ED on a a/8 x \/8 cluster for pdn/pda = 
—0.2. (a) Results corresponding to several U's and J/[/ = 
1/8. (b) Results varying J, at fixed U = 1. 
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FIG. 11: (Color online) Spin structure factor S(k) for (a) 
half filling and pdn/pda = 0.1, U = 0.5, J = C//8, and (b) 
pdn/pda = -0.2, [/ = 0.5, J = i7/8 at half filling and with 
two additional electrons. These are ED results on a \/8 x \/8 
cluster. 



pcctcd, rather than genuine long-range order). Increas- 
ing the Hund's coupling J at a fixed U, see Fig. |lO(b)[ 
produces a similar effect, because it leads to larger local- 
ized moments, allowing for a stronger overall collective 
spin ordering. However, for positive pdn, the diagonal 
hopping ^3 is no longer strong enough to drive the (tt, 0) 
order, and the spin structure factor peaks at (tt, tt) in- 
stead (Fig. 11(a) I . Figure 11(b) shows the spin structure 
factor for pdn /pda = —0.2, U = 0.5, and J = U/8 when 
two more electrons arc added to half filling. The (0, n) 



order is weakened in the doped system. 

In our previous effort,— we investigated the pairing 
symmetry for two added electrons in the region of hop- 
pings —0.5 < pdn < —0.2, varying U, and for the special 
case J/U = 1/4. The spin of the state with two ad- 
ditional electrons can be determined by comparing the 
ground state energy for a total z component of the spin 
Sz = to Sz = 1- If these two energies are degenerate, 
the state is a triplet (it was also tested that the ground 
state of 5^ = 2 is not degenerate with 3^ = and 1, 
thus excluding higher spin states). The Hubbard repul- 
sion U was found to drive the spin of the two electrons 
added to the half-filled system from triplet at small U 
to singlet at larger U, for the pdn^s investigated. The 
critical Uc needed for the transition was found to be the 
lowest at pdn « —0.2. This value of pdn moreover leads 
to a Fermi surface with hole and electron pockets similar 
to that obtained with band structure after folding [see 
Ref. I?! and Fig. |9(b)] . For large \pdn /pda\, on the other 
hand, the Fermi surface has far larger electron pockets 
around (0, tt) and (tt, 0) than th ose fo und in band calcu- 
lations or experiments, see Fig. 9(a) Consequently, we 
will mainly focus on pdn/pda = — 0.2 below. 

In Fig.[T^ the regions where singlet and triplet pairing 
dominate, depending on U and J, are shown. (The no- 
tation "singlet 9" and "singlet 2" refer spin-singlet states 
with and Ai^ symmetry, as discussed in more de- 
tail in the "Pairing symmetry" section below as well as 
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FIG. 12: (Color online) Dominant pairing tendencies of the 
ground state corresponding to two more electrons than half- 
filling, at pdTT=-0.2. The notation "singlet 9" and "singlet 2" 
is explained in the text. 



in App. [B1) The trends observed for J/U = 1/4 re- 
main stable for other realistic values of J. Additionally, 
qualitatively we have observed that increasing the Hund 
coupling J promotes a robust triplet pairing. Due to 
the relation used between U, U' , and J, two electrons 
on the same site, but in different orbitals, no longer feel 
any Coulomb repulsion for the maximal J = 2U/5. As 
a consequence, values oi J/U > 0.4 are here considered 
unphysical. 



2. Pairing symmetry 

To investigate the symmetry under rotations of the 
pairing states, the half-filled ground state is compared 
to states with two additional electrons. The symmetry 
sector of the half-filled ground state must be contrasted 
with the symmetry of the doped state, and the symme- 
try operation leading from one to the other gives us an 
indication for the pairing symmetry. This method was 
very successful in establishing the c?-wave character of the 
pairing in the t-J model for the cuprates.— In addition, 
we also added a pair of electrons with a well-defined sym- 
metry under rotations to the half-filled ground state, and 
calculate the overlap between the resulting state and the 
ground state obtained for half filling plus two electrons. 

From this analysis, we found that the dominant pairing 
operator for spin-singlet pairs is inter-orbital and given 
byZa 

At _ I V (d^ -d^ d^ ) 



(23) 

where i = 1, . . . , iVgitcs denotes the lattice site, jl ~ x,y 
the unit vector connecting NN sites, and a = x,y the 
xz and yz orbitals, respectively. This operator trans- 



FIG. 13: (Color online) Overlap between the ground state for 
half filling plus two electrons \(j>2), and the state obtained by 
applying the pairing operators and Ag to the undoped 
ground state \<j)o}- At small to intermediate U, the inter- 
orbital singlet pairing operator AJ with B2g symmetry {i = 9) 
dominates, see Eq. (|23|) . For large U, the intra-orbital singlet 



Aj {i~2) dominates. Results were obtained using pdn/pda — 
—0.2, and the ED technique on \/8 x \/8 site clusters. 



forms as S2g, and it is #9 in the detailed list provided 
in Ref. [H. In addition to the pairing between near- 
est neighbor sites, we also find a small overlap for the 
corresponding onsite pairing #8 (reaching at most 
10% of the intersite overlap) and some overlap for the 
NNN i?2g pairing. In contrast to the small onsite contri- 
bution, the NNN pairing is sizable, but its exact strength 
compared to NN pairing is difficult to ascertain with the 
small cluster used. 

The only other singlet pairing for which we have found 
a substantial overlap is #2 in the above mentioned list, 
although, as discussed below, its region of stability at 
large U does not have the correct properties expected for 
the FeAs new superconductors. This operator is intra- 
orbital, has Aig symmetry, and it is given by 



1 



). (24) 



Applying the pairing operators A- to the half-filled 
ground state |^o)i we find that the resulting vector 
a] |(/)o) has a very small norm < 0.15 for pairings i = #3, 
#4, #5 and #6 in the list of Ref. [s^, while it reaches 
w 0.6 - 0.8 (depending on U and J) for #1, #2, #7, 
#8 and #9. Then, at least qualitatively, we conclude 
that only the latter pairs can be created easily in the 
half-filled ground state. To provide more quantitative in- 
formation, we then calculate the overlap between A||(/)o) 
and the ground state found for half filling plus two ad- 
ditional electrons \(p2)- We only find substantial over- 
laps for the operators #9 (i52g) and #2 (^ig) given in 
Eqs. ([221) and ([M]) . while Ay|0o) is always orthogonal 
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FIG. 14: (Color online) Spin structure factor S{k) for hopping 
parameters ti = —1.0, t2 = 1.3, ts = = —0.85 (in eV 
units) "2^ (a) results for several values of U and J/U = 1/8; 
(b) results for two values of J, with U = 2.8 eV fixed, (a) 
and (b) were obtained at half filling using ED on X 78 
clusters. 
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FIG. 15: (Color online) Spin structure factor S{k) for half 
filling and half filling plus two electrons, using the hopping 
parameters ti = —1.0, t2 = 1.3, ta = t4 = —0.85,— and 
[/ = 2.8, J = 0.1 (in eV units). These are ED results for 
VS X \/8 clusters. 



to the two-electron ground state \(f>2), at least for the 
range of parameters investigated. As it can be observed 
in Fig.[T21 the i?2g pairing Eq. (P5|) occurs at small to in- 
termediate Coulomb repulsion U, which is the expected 
suitable regime to describe non-insulating materials with 
bad metallic properties. Only for large U ^ 2.8 eV, where 
a hard gap in the density of states indicates insulating be- 
haviorf2& we do find the pairing Eq. with Aig sym- 
metry. In this regime we find some admixture of the 
corresponding onsite pairing (#1) and the longer-range 
j4ig NNN pairing which, as shown in App.lBl corresponds 
to the much discussed .s± pairing statej ^^i^^i^^i^^ 

Figure [12] shows more explicitly the regions of dom- 
inance of the states #9 (i32g) and #2 (Aig) in the U 
vs. J/U plane. Thus, the i?2g symmetric operator seems 
to be the most realistic in the regime pdir/pda w —0.2 
and 0.5 < J7 < 1, which is the appropriate region of 
parameters to qualitatively describe the new FeAs-based 
superconductors 
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FIG. 16: (Color online) Dominant pairing tendencies of the 
ground state for two more electrons than half-filling, with 
hoppings fitted to band-structure results?^ The "singlet 9" 
and "singlet 2" notation is the same as for the case of the 
Slater-Koster hoppings. 



C. Results with hopping parameters fitted 
to band-structure calculations 

1. Results at nonzero J 

We have also investigated the two-orbital model using 
hopping parameters obtained from a fit to band-structure 
calculation resultsi^ It is interesting to observe that this 
set of parameters also leads to (0, tt) / (tt, 0) antiferromag- 
nctic order at half filling, see Fig. [Ml which is again en- 
hanced by increasing U at fixed J, or increasing J at 
fixed U. As for SK hoppings. Fig. [T5| shows that the 
magnetic order is only slightly reduced by the doping 
with two electrons. In Fig. I16[ we report the spin of 
the state with two electrons added to the half-filled state 
and find qualitatively similar behavior as with the SK 
approach: U promotes singlet pairing, in the previously 



discussed "#9" (i?2g) and "#2" (^ig) channels, and J 
favors triplet pairing. 

We have performed an analogous analysis of the elec- 
tron pairing as in the previous subsection, and again find 
the inter-orbital singlet operator Eq. ([^5)) to dominate at 
intermediate U and J. Table |TTT| gives the pairing ampli- 
tudes for operators Eqs. (|^5)) and (|24p . for several ([/, J) 
parameter sets. As before, the Coulombic parameters 
were chosen to give a spin-singlet state for two electrons 
added to half filling, but are expected to be small enough 
to remain in the metallic regime. As for SK hoppings, 
the pairing symmetry for these singlet states is i?2g, i-C, 
the inter-orbital singlet Eq. (|^5)) dominates here as well. 

Inter-orbital pairing is favored over intra-orbital pair- 
ing by the kinetic energy, because the inter-orbital 
Coulomb repulsion U' is weaker than the repulsion U 
within the orbitals. To test this assumption, we analyze 
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TABLE III: Overlap between the ground state 102) for two 
electrons added to half filling and the states A||0o) that are 
obtained by applying pairing operators Eqs. (|23|l and (|24|l to 
the ground state at half filling. Data are for a kinetic energy 
operator obtained from band-structure fitting with ti — —1.0, 
t2 = 1.3, and ts = t4 = -0.85 (in eV units).— 
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TABLE IV: Overlap {(f>2\Al\(f>o) as in Tab.|TITl but with J = 0. 
The coupling U' is different from U, and the hoppings were 
ti = —1.0, t2 = 1.3, and ts = t4 = —0.85 (in eV units). 
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pairing amplitudes for the special case J = but instead 
of using U — U' , we use U' < U, i.e., we deviate from 
the relation [/' = [/ — 2 J. The resulting pairing symme- 
try is still the same B2g inter-orbital singlet Eq. as 
deduced from the amplitudes on Tab. |TVl 



2. Results at J =0 and U = U' 

Finally, let us discuss the special case J — 0, U — U' 
for U = 2.8 (in eV units), which are the couplings used 
in Ref. [zl. In contrast to J > 0.05 and U' < 2.6, we here 
find that the Ai^ intra-orbital singlet Eq. ([M]) has the 
lowest energy, i.e. the same pairing as observed at large 
U for both the fitted and the SK hoppings. However, at 
L/ = [/' = 2.8 and J = 0, the second lowest eigenstate 
is almost degenerate with the ground state, and it gives 
an overlap with Eq. (pS)) . i.e. pairing #9. The third 
state at these couplings corresponds to the intra-orbital 
pairing #7 (i^ig). Table IVl contains the explicit numbers 
for energies and overlaps. The near degeneracy of these 
states does not allow us to reach a clear conclusion for 
the J = and U = U' case, which appears to be singular, 
since small modifications away from our results, such as 
increasing lattice sizes, may change the relative order of 
the competing states. 



TABLE V: Energy of the lowest eigenstates for two electrons 
and J — 0, U = U' — 2.8, and pairing operators giving the 
largest overlap when applied to the half-filled state. Hoppings 
are those from band-structure fitting. 



energy 


pairing 




1 -8.45322 


#2 


0.66 


2 -8.45150 


#9 


0.66 


3 -8.4132 


#7 


0.67 


TABLE VI: Four-spin correlations (v^^lvx) and {vy\vx) for two 
parameter sets using the full two-orbital model at the values 
of U and J indicated, on an 8-site cluster. For comparison, 
results using a perfect Ising-like (0, vr) + (tt, 0) state on 8 sites 
are also shown. 


Ising 


[7 = 4 [7 = 1 


U = 0.5 




J = 1 J = 0.25 


J = 0.0625 


{Vy\v^) -1 


-0.64 -0.1684 


-0.07776 


{v^\v^) 0.0625 


0.02735 0.02856 


0.0225 



3. Competing magnetic states at half-filling 

The magnetic phase diagram for the present model was 
also studied using a mean-field approximation!^ At half 
filling, it was claimed that the Coulomb repulsion rather 
than stabilizing a state with (tt, 0)/(0, tt) spin-stripe or- 
der, induces an "orthomagnetic" (OM) ordering where 
NN spins are at right angles. The spin structure factor 
for this phase is still peaked at (0, tt) and (tt, 0). We have 
tried to address this issue with the ED technique on the 
8-site cluster. Unfortunately, several observables are ex- 
pected to give similar results for the two states. Apart 
from having similar spin structure factors, diagonal NNN 
spin-spin correlations are also negative in both states. 
Moreover, even the expectation value for NN correlations 
vanishes in both cases: in the OM state because the NN 
spins are at 90°, while in the stripe state the cancella- 
tion occurs because the small cluster ground state |(/)o) 
contains both (0,7r) and (tt, 0) configurations with equal 
weight and NN correlations {4>o\vx) and (0oky) average 
out, where 

\v.)^^J2^i-S,+M , (25) 

jVsitcs . 

1 

\vy) = ^J2^i-Si+M ■ (26) 

-'Vsitcs 

1 

Our numerical ED results indeed give very small negative 
values for the NN spin correlation. However, we expect 
the two states to lead to different results for {vx\vx) and 
{vy\vx)- {vx\vx) is expected to be positive and {vy\vx) 
negative in the spin-striped phase, while both should be 
zero or very small in the OM phase. Table IVTl shows the 
results for strong (U = 4), intermediate {U = 1), and 
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weak {U = 0.5) on-site Hubbard repulsion. For [7 = 4, 
we clearly find {vy\vx) < and {vx\vx) > 0. While 
these numbers become weaker for the less spin ordered 
states at smaller U, they are still consistent with the 
(0, 7r)/(7r, 0) ordering. For comparison, we also include 
results obtained for a linear combination of states with 
perfect Ising-like (0, tt) and (tt, 0) order in the z-direction. 
These results are useful to judge the expected order-of- 
magnitude values for the correlations investigated. The 
numbers on Table VI show that the numerical results 
for the two-orbital model are compatible with those of 
the Ising spin-stripe state, particularly considering that 
quantum fluctuations will reduce the spin correlations 
of such a state. As a consequence, our present inves- 
tigations favor the spin-stripe magnetic state, although 
further work is needed to fully confirm these conclusions. 



V. DISCUSSION OF NODAL STRUCTURE IN 
THE MEAN-FIELD APPROXIMATION 

In this section, the results of a pairing mean-field anal- 
ysis of the two-orbital Hamiltonian will be discussed. The 
numerical results of the previous sections and experimen- 
tal data will be used to guide this mean-field approxima- 
tion. Experiments indicate that the pairs in the Fe-based 
superconductors are spin singlets^i^i^S and our previ- 
ous numerical results did provide a dominant spin-singlet 
pairing operator, as discussed in the previous section. It 
is also important to notice that the pairing operator that 
we obtained mixes different orbitals. A numerical study 
of the orbital composition of the bands that determine 
the FS in our two-orbital model indicates that the bands 
that constitute the pockets are an admixture of xz and 
yz orbitals. Thus, it is not surprising that the dominant 
pairing operator creates pairs made of electrons in dif- 
ferent orbitals. One of the main goals of the analysis 
discussed below will be to find out the nodal structure of 
the mean-field Hamiltonian. 



A. Location of the Nodes 

1. Reminder of one-band model results. 

For the simple case of c?-wave superconductivity in a 
single-orbital model, characterized by the dispersion re- 
lation ^(k) = — 2t(cos kx + cos ky) — At' cos cos ky — /i, 
the gap function is given by = A (cos kx — cos ky). In 
this case, the mean- field Hamiltonian reduces to a 2x2 
matrix linking k with -k which is simply given by 



values that solve the equation det{HMF) = 0, i.e., ^(k)^ 



(27) 



To obtain the position of the nodes in the gap, we merely 
need to find the values of k^ and ky where the eigenvalues 
of the matrix Eq. (P7)) are zero. These are the same 



0, which is satisfied only if each term vanishes 



independently. This occurs at the points where the non- 
interacting Fermi surface described by ^ (k) = intersects 
the diagonal lines along which A^ = 0, i.e., k^ = ky and 
kx = ~ky. This procedure establishes the well-known 
location of the four d-wavc nodes of a single-band model. 



2. Nodes in a two-orbital model 

For a system with two orbitals, we will proceed in an 
analogous manner as for one orbital. The MF Hamilto- 
nian matrix in the basis (dj^. x p^ky p '^-k,x,i, <^-k,i/,i) is 
now given by the 4x4 matrix 



Hmf — 



S,xy Cyy 

Afe 
VAfc 







Afc 


~^xy 



(28) 



The matrix elements can be obtained by Fourier trans- 
forming the tight-binding Hamiltonian i?TB given in 
Eq. HH). We obtain: 



^xx — 

^xy = 

and 



—2t2 cos kx — 2ti cos ky — At^ cos kx cos ky — /i, 
—2ti cos kx — 2^2 cos ky — At^ cos kx cos ky — 
— At 4 sin kx sin ky, (29) 



T^(cos kx + cos ky), 



(30) 



where V is the strength of the pairing interaction. 

Notice that we can also work in the basis in which 
-ffxB is diagonal. In this basis, which is expanded 

by (4,i,T'4,2,T'C-fe.2a'C-fe,i,i), ff] 

U^^HmfU given by: 



MF becomes iJ^F = 



/ £1 Vb Va\ 

e2 ~Va Vb 

Vb -Va -62 

\Va Vb -ei / 



(31) 



where Va and Vb are given by: 

Va = 2uvAk 



Vb = {v^~u^)Ak, 



(32) 



(33) 



and u and v are the elements of the change of basis matrix 
U given by 



U 



( u V \ 

w -u 

u 

V -u V ] 



(34) 
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FIG. 17: (Color online) (a) Energy vs. momentum for the 
non-interacting tight-binding Hamiltonian Eq. (|14|l using the 
hoppings from band-structure calculations,— ti = —l,t2 — 
1.3, and ts — t4 = —0-85 (in eV units). These non- interacting 
results are shown as continuous lines. Also shown are the ad- 
ditional Bogoliubov bands produced by the pairing interaction 
considered in this work (dashed lines). The circles (boxes) in- 
dicate the bands that contribute electrons to the intraband 
(interband) pairs at different locations close to the FS (for a 
discussion see text), (b) Fermi surface for the corresponding 
non-interacting half-filled system. 



with U ^ ~ . Remember that Va, Vb, u, and v are 
all functions of the momentum fc, and + v'^ = 1. 

Now let us discuss the physical meaning of Va and Vb ■ 
According to Eq. (|3T|) . Va is the intraband pairing for 
band 1, i.e., the band with the highest energy (electron 
band), while the intraband pairing for band 2 (hole band) 
is —Va- Thus, there is a relative phase tt between the two 
intraband order parameters. In the standard BCS studies 
for multiband models, it is expected that pairs are formed 
by electrons in the same band.— In the two orbital model, 
as just discussed, we found intraband pairing but we also 
obtain interband pairing with strength Vb ■ The possibil- 
ity of interband pairing has been considered previously 
in several contexts: (i) Possibility of T^, (ii) high Tc 
cupratesr^ and (iii) heavy fermion systems^Si where it 
was shown that if two bands are very close to each other 
in the vicinity of the Fermi level, interband pairing can 
occur. The weaker the pairing potential the closer to the 
Fermi surface the two bands have to be. When long range 
pairing develops the Brillouin zone gets folded and, as a 
result, the total number of bands doubles. In this rep- 
resentation, which arises by diagonalizing Eq. (|28|) with 
Afe = 0, the two-orbital model has the dispersion shown 
in Fig. [TTI where each band (panel (a)) and FS (panel (b)) 
is represented with a different color and the folded (un- 
folded) portions with dashed (continuous) lines. It can be 
seen from panel (a) that at the Fermi level there is only 
intraband crossing indicated by green circles. However, 
there is also interband crossing, indicated by the orange 
boxes, above the Fermi energy. The previous numerical 



results appear to indicate that the effective coupling is 
sufficiently strong as to produce interband pairing. In 
Fig. [18] it can be seen that even a small V opens a gap 
between the two bands that cross away from the FS. In 
multiorbital models the opening of these gaps can lower 
the overall energy.iiS^ Within the standard BCS approach 
this result may appear counterintuitive and it could be 
an artifact of the two-orbital model or of the small sys- 
tem size that we have considered. However, there are 
clear indications that most of the FS in the five-orbital 
model have character d^z and dy^^^ and, in such a case, 
the only possible pairing operators that respect the sym- 
metry of the FeAs planes are those that have been con- 
sidered in our calculations and othersi^ In fact, none of 
the 16 pairing operators that combine these two orbitals 
leads to a purely intraband pairing interaction. This, of 
course, could be an indication that other orbitals have 
to participate in the model but we believe that it is still 
instructive to consider the nodal structure that results 
from the pairing operator that was favored numerically 
within the two-orbital model and attempt to compare the 
results with experiments. 

The existence and position of nodes in the resulting 
mean-field band structure can be found by requesting 
that det(i?^p) = 0, as in the one-orbital case. From 
Eq. (pij). we obtain the following equation: 

VliVl + 2Vi +el + 4) + (£162 + Vif = 0. (35) 

This equation is satisfied in two possible ways: 

(1) First, a solution can be found if Va = Vb = 0, 
and ei = or £2 = 0. Thus, this condition for nodes is 
satisfied if the lines where cos + cos ky = 0, namely 
the lines where Va and Vb vanish, intersect any of the 
non-interacting Fermi surfaces determined by the points 
where ei = or £2 = 0. It is clear that the line 
cos kx + cos ky = intersects each of the four electron- 
pocket Fermi surfaces in two points per pocket (see 
Fig. [T7b and Fig. le). This means that nodes will ap- 
pear only in the electron pockets, not in the hole pock- 
ets. These are the nodes that arise from a simple ex- 
trapolation of the reasoning used to find nodes in the 
one-orbital model, namely by finding the intersections of 
the non-interacting Fermi surface with the trigonometric 
function, in this case cos k^ + cos ky , contained in the Afc 
gap function. The position of the nodes in the electron 
pockets upon folding of the Brillouin zone (see Fig. Id 
and Fig. [TTb) is at the points in fc-space where the two 
electron pockets intersect each other. 

Notice that the existence of these nodes does not de- 
pend on the value of V. They will always be present as it 
can be seen in Fig. [T8l where the nodes along the X — Y 
direction appear in all the panels varying V. 

(2) However, the two-orbital nature of this problem 
leads to the possibility of additional nodes in unexpected 
locations. This can be understood by realizing that 
Eq. ((55l) can also be satisfied if Va = and 

Vb' = -£i£2. (36) 
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FIG. 18: (Color online) (a) Energy vs. momentum for the 
mean-field Hamiltonian using ti = —l,t2 = 1.3, and — 
t4 — —0.85 (in eV units).— The Bogoliubov bands produced 
by the pairing interaction considered in this work are also 
shown, with equal intensity. The four panels correspond to 
four different values of the pairing attraction: (a) V = 0, (b) 
V = 0.5, (c) V = l, and (d) V = 8 



According to the expression for Va in Eq. ([32|) . and as- 
suming that Afe is non-zero (if it is zero we recover the 
nodes already described in (1) above), then the condition 
Va = is satisfied ii u — or v — 0. Due to the nor- 
malization + v"^ ~ 1, when u ~ then it must occur 
that V — 1, and vice versa. Introducing these values of 
u and V in Eq. ([34|), it can be shown that the condition 
that i?MF " U^^HmfU is diagonal is satisfied only if 
£,xy = 0- According to Eqs. (|29p . for S^^y to vanish it is 
necessary to have = or tt, or ky = or tt. Then, 
new nodes could be expected along these horizontal or 
vertical lines in the Brillouin zone. Since the product of 
the two energies £162 has to be negative (i.e. the energies 
cannot vanish, otherwise we recover (1)), the nodes, if 
they exist, will appear in between the hole and electron 
pockets at locations in fc-space that do not belong to the 
original tight-binding Fermi surface. To understand this 
interesting result, consider the example of k^ = 0- For 
this special case we obtain, 

Vi = V\l+cosky)^, 
ei = — 2f2 — 2ii cos fcy — 4i3 cos fcy — ^, 



£2 



-2ti — 2t2 cos ky — Ats cos ky 



(37) 



Replacing Eqs. ([37]) in Eq. ([36]) . a quadratic equation is 
obtained that allows us to find the values of cos ky where 
nodes should appear. Depending on the specific values of 
V, the hopping amplitudes, and ^, the quadratic equa- 
tion can have two solutions (meaning that two nodes ap- 
pear along the k^ = axis between the hole and electron 
pockets), or one solution (meaning just one node), or no 
solution at all (indicating no extra nodes). Thus, once 



the folding and rotation of the FBZ is performed, nodes 
can appear along the diagonals of the BZ in Fig. Ic for 
particular values of the parameter in the model.^^ 
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FIG. 19: (Color online) Examples showing the positions of 
the nodes for the two-orbital model, using the standard pair- 
ing mean-field approximation. With solid lines are shown 
the hole pockets (red) and electrons pockets (green), namely 
the Fermi surface in the non-interacting limit V" = 0. The 
black crosses are the nodes. As example, we use the case of 
the hopping amplitudes derived from fits with band-structure 
calculations.— (a) corresponds to the "weak" pairing regime, 
V = 0.5, showing the existence of many nodes, both at the 
electron pockets as well as near the hole pockets, (b) is the 
example of V = 2.0 where the number of nodes has been re- 
duced compared with (a), (c) is the "intermediate" V = 3.5 
regime {V is still substantially smaller than the bandwidth) 
that provides the smaller number of nodes (a total of 8), all 
being located at the non-interacting electron pockets, (d) is 
the folded Fermi surface corresponding to case (c) for better 
comparison with experiments and band calculations. 



A variety of examples obtained numerically illustrate 
this nontrivial nodal structure, as shown in Fig. [191 (a) at 
weak V, several nodes are found either at or close to both 
the hole and electron pockets. In view of recent photoe- 
mission experiments reporting the absence of nodes at 
the hole pockets (see discussion below), this regime is 
unlikely to be realized experimentally, (b) is obtained 
increasing V: in this case the number of nodes has de- 
creased. In addition to those coming from solution (1) in 
the previous discussion, all at the electron pockets, still 
solution (2) provides some nodes at the boundaries of the 
Brillouin zone in this regime, (c) is the most canonical re- 
sult, obtained at intermediate V, with the nodes only ap- 
pearing in the electron pockets where cos k^ + cos ky = 
intersects the original Fermi surface. Both in (b) and (c) 
there are no nodes in the T centered hole pocket even for 
this i?2g state, (d) provides the results of (c) but in the 
folded zone for comparison with experiments. 

Band dispersions obtained with mean field results are 
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also shown in Fig. [THl We observe how the nodes along 
the r — X and X — M directions get closer to each other 
as V increases from to 0.5 and to 1; and how they have 
disappeared for V=%. It is also interesting to see how 
the crossing of different bands, indicated by the orange 
squared boxes in Fig. [17] (a), is replaced by a gap as soon 
as V is finite (sec panel (b) in Fig. [T5)) which appears to 
be the effect of the interband interaction. 

For completeness, in Fig. [201 we provide the nodal 
structure for the case of hoppings obtained from the 
Slater Koster approximation, that gives large pockets in 
the band-structure calculations. In the weak coupling 
case, (a), once again several nodes are obtained. This 
regime appears unrealistic. Increasing V , panel (b) shows 
that the nodes only remain in the electron pockets, as 
found before in Figs. [T9| (c,d). 



(a) (b) 
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tron pockets were also studied. However, in that effort 
the 122 materials for which the FS depends on kz were 
analyzed, and only a few cuts in momentum-space were 
investigated. Recent ARPES experiments that suggest a 
short Cooper-pair sizoi24 would suggest that the regime 
of large V in our study of the nodal structure is the most 
realistic, thus clearly locating all the nodes in the electron 
pockets. Then, a more detailed ARPES analysis would 
be needed to fully conclude that there are no nodes in this 
system, particularly in the electron pockets f^S^ Other re- 
cent ARPES experiments have shown a variety of inter- 
esting aspects, such as substantial differences with band 
structure calculations that also need to be incorpo- 
rated in future theoretical studies. 

The information provided in this manuscript for 
the actual location of the nodes for the i32g state 
will help to guide future ARPES experiments. In 
view of the several other experimental investiga- 
tions that have reported nodes in the Fe-based 
superconductors^^'^^'^^^^i^^ we believe 
this issue is still open and needs further research to ar- 
rive to a final conclusion. If future experimental work 
clearly proves that there are no nodes in the new Fe- 
based superconductors, not only in the hole pockets but 
more importantly in the electron pockets, then it will be 
concluded that the two-orbital model used here will not 
be sufficient to properly describe this family of materials, 
and more orbitals will be neededj^Si 



FIG. 20: (Color online) Nodal structure for SK hoppings with 
(pdn) — —0.2. (a) is in the weak coupling regime, V — 0.05, 
and shows many nodes as in Figs. 1191 fa. bl. (b) is the interme- 
diate coupling regime, V — 0.5, where the number of nodes is 
the minimal (8), similarly as in Figs. I19l (c.d). 



In addition to the analytic discussion, we have also 
searched numerically, using a large lattice 200x200 in 
fc-space, for the zero eigenvalues of the original matrix 
Eq. pS)) . These numerical results are in excellent agree- 
ment with the analytic discussion, thus showing that 
the nodal structure of the two-orbital mean-field pairing 
Hamiltonian has been properly obtained. 



3. Comparison with ARPES experiments 

How do these theoretical calculations based on the two- 
orbital model compare with angle-resolved photoemission 
experiments for the Fe pnictides? In Ref. |4^, ARPES re- 
sults were presented with the focus of the effort on the 
hole pockets at F. It was concluded that nodes were not 
observed in those hole pockets. This result is compatible 
with our state since nodes do not appear on the hole 
pockets, but instead on the electron pockets, at least at 
intermediate values of the attraction V. In Ref. a 
similar conclusion was reached but in this case the elec- 



VI. CONCLUSIONS 

In this manuscript, we have studied some properties of 
a two-orbital approach for the new Fe-based supercon- 
ductors. It is important to find out the minimal model 
capable of reproducing the basic physics of these materi- 
als. By studying a relatively simple model, considerable 
insight could be reached on the inner mechanisms that 
cause magnetism and superconductivity in these com- 
pounds. While models with more than two orbitals would 
certainly be more accurate, the difficulty in extracting 
reliable numbers from the models grows fast with the 
number of orbitals. 

Here we have shown that the magnetic properties of the 
undoped parent compound are properly reproduced by 
a simple two-orbital model: spin stripes are obtained in 
agreement with neutron scattering experiments. Regard- 
ing electron doping, here we follow the same approach as 
for the cuprates: it is expected that the pairing channel 
will be unveiled by simply studying the symmetry prop- 
erties of the state of two electrons added to the half-filled 
ground state. This approach worked for the models for 
Cu-oxides superconductors, leading to the d-wave state 
prediction. Within this assumption, the spin-singlet pair- 
ing state that dominates in the phase diagram at realistic 
values of the Hubbard repulsion U is found to transform 
according to the i?2g representation of the lattice sym- 
metry group. At large Coulomb repulsion U, too large 
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to describe the metallic state of the undoped compound, 
we found that the relative symmetry of the undoped and 
electron doped ground states is the same and, thus, they 
are connected by a pairing operator that transforms ac- 
cording to Aig. We showed that the NNN pairing oper- 
ator with this symmetry is the "s±" state and that this 
state, according to our numerical calculations, prevails 
only in an unphysical regime of parameters. On the other 
hand, for a robust electron-electron effective attraction to 
form Cooper pairs, assumption compatible with the con- 
clusions of recent ARPES experiments, the pairing 
state found for realistic U has nodes only in the electron 
pockets. All our main conclusions do not depend qualita- 
tively on the set of hopping amplitudes used: our results 
appear to be representative of the two-orbital framework 
in general and not merely of a particular model with par- 
ticular couplings. Thus, a conclusion of our study is that 
more refined ARPES experiments in the superconduct- 
ing state are needed to analyze the possible existence of 
nodes in the electron pockets. These future experiments 
will provide crucial information to guide the theoretical 
search for the minimal model that captures the physics 
of the Fe pnictides. 
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APPENDIX A: EFFECTIVE INTERACTION 
THAT GENERATES THE Bag PAIRING 
OPERATOR 

The mean-field superconducting state discussed be- 
fore could originate from an effective attractive density- 
density interaction dynamically generated in the original 
Hamiltonian, or induced by particular phononic modes if 
an electron-phonon coupling is incorporated. The form 
of this attraction is: 

^fattr = 51 '^i,a,ani+fj.-a,-a, (Al) 

and below we prove that indeed it generates the correct 
pairing term. It is well-known that a similar nearest- 
neighbor density-density attraction of the form -Vnini^fj^ 
leads to d-wave superconductivity in a mean-field treat- 
ment of the one-band repulsive U Hubbard model^^ and 
here we merely generalize this concept to two orbitals. 

Let us discuss the mean- field treatment of -ffatti - In 
momentum space, the Fourier transformed of this effec- 



tive attraction is: 

k,k'.a,cr 

(A2) 

where we have requested that the pairing occurs be- 
tween electrons with opposite momentum (thus, we have 
dropped a third sum over all wavevectors), in different 
orbitals, and with opposite spin, as required by the dom- 
inant singlet pairing operator obtained from the numer- 
ical simulations. The potential is given by 

Vk,k' = -2V[coa{k'^ ~ k^) + cos{k'y - ky)] 
= F(k'-k)=^Kr/,(k)r;,(k'), 

i 

where 77i(k) are the irreducible representations of the 
group . Since the ED numerical results indicate that 
the pairing operator is proportional to (cos/c^. -I- cos/cy), 
which corresponds to the irreducible representation Aig, 
we will focus on that particular term in the expansion of 
the full potential Vk,k' Thus, we will consider Eq. 
but using the long-range separable potential 

Vic.k' = V*{coskx + cos fcy)(cos fc^ + cos k'y), (A3) 

instead of the full short-range potential. 

We will treat the four-fermion term in i?attr within 
the usual mean-field approximation^^iS where some pairs 
of fermionic operators are replaced by numbers, such as 
(^k a)' found self-consistently. Then 

Hmf = Htb + ^ (^k,k'(&k_Q)<^-k',-Q,irfk',Q,T + 

k,k' ,a 

Vk,k'(fek',a)<„,T^Lk,-.a)- E ^k.k'(<a)(&k',.). 

k,k' ,o 

(A4) 

Defining 

A(k) = -^yk,k'(fek',a), 

k'jQ 

At(k) = -^Vk,k'(6L,J, (A5) 

k',a 

and using the separability of the potential Eq. (|A3p we 
obtain 

Hmf = HTB + J2iA^i^)d^K (A6) 

k,Q 

A(k)<„,^dlk,-aa) - E ^k,k'(&L)(6k',a), (A7) 

k,k',a 

which leads to the same self-consistent equations as in 
Section V by setting 

At(k) = A(k) = V* A{coa k^+ COS ky), (A8) 

and V ~ V*A, where A should be obtained by solv- 
ing the gap equation that is obtained from minimizing 
the energy of the mean-field Hamiltonian with respect to 
A(k). 
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APPENDIX B: s± PAIRING INVOLVING dx 



In Section IV. B. 2 we showed that our numerical sim- 
ulations favored a spin singlet interorbital pairing state 
with symmetry in the physical regime of parameters 
of the two orbital model, while a pairing state with sym- 
metry Aig prevailed only in the unphysical strong cou- 
pling regime and at the singular point U = U' , J — 0. 
In this appendix we will discuss in more detail the pair- 
ing operators with Aig symmetry in the context of the 
two orbital model, and we will show that a s± pairing 
operato r ^^i^^i^^i^^ involving only d^z and dy^ electrons be- 
longs to this group. 

The extensive literature on the s± pairing 
state^i^i^'^ indicates that this state does not 
have nodes on the Fermi surface and that 



Al(k) = -A2(k + q) 



(Bl) 



where 1(2) denotes the electron (hole) Fermi surface and 
q= (7r,0) or (0,7r),^ 

From the classification of possible pairing states for 
the dxz and dyz orbitals provided in Ref. |89| we realize 
that there exist the following four nodeless pairing oper- 
ators: (i) pairing state #1 which produces on-site intra- 
band pairs which are even under orbital exchange, spin 
singlets, and transforms according to the Aig irreducible 
representation of D4h. Following the steps of Section V.A 
for this pairing state we obtain: 



MF 



^ £,xx S,xy Aq 

^xy ^yy Aq 

Ao ^^xx ^£,xy 

V Ao -ixy ~^yy J 



(B2) 



where Aq is a constant independent of momentum. In the 
base in which the tight binding Hamiltonian is diagonal 
we obtain: 



H' 



MF 



/ ei Ao \ 

62 Ao 

Ao -£2 

V Ao -ei / 



(B3) 



This leads to intraband pairing Aq = V which is momen- 
tum independent and equal for the two bands. This does 
not correspond to the s± pairing since it does not satisfy 
Eq. ((BT|) ,^ 

Now let us consider nearest-neighbor pairing. We find 
that the only nodeless pairing operators involving elec- 
trons in nearest-neighbor sites also have symmetry Aig 
and result from a (ii) symmetric (or (iii) antisymmetric) 
combination of pairings #2 and ,^3. The symmetric (an- 
tisymmetric) combination corresponds to pairing of the 
dxz electrons along the x (y) direction while the dyz pairs 
along the y (x) direction. Following the steps of Section 
V.A we obtain: 



MF 



^ £,xx ixy Ai \ 

S,xy S,yy A2 

Al ~£,xx —£,xy 

\ A2 -S,xy -iyy J 



(B4) 



where Ai = T^cosfc^, {V cos ky) and A2 = V cos ky 
{V coskx) for the symmetric (antisymmetric) combina- 
tion. In the base in which the tight binding Hamiltonian 
is diagonal we obtain: 



^MF — 



/ ei V12 Vi \ 

62 V2 V12 

V12 V2 -62 

V Vi V12 -61 / 



(B5) 



where Vi, V2 and V12 are given by: 

Vi = u^Ai +w^A2, 

V2 = w^Ai 4-u^A2, 

V12 = uw(Ai - A2). 



(B6) 
(B7) 
(B8) 



Thus, this leads to intraband interactions Vi and V2 
which, according to our numerical checks, satisfy l/i(k) = 
— Vi(k-l-q) as expected for the s± pairing^ but there 
is interband pairing given by V12 which is considered un- 
physical by many authors. ^^^iSiiS^iS^ 

(iv) Finally, we can also focus on pairs of electrons 
along the diagonals of the square lattice formed by the 
Fe ions. Following the notation of Ref. [s^ the corre- 
sponding basis function is cos kx cos ky that transforms 
according to Aig. This basis provides a pairing opera- 
tor with a full gap and which is a spin singlet. It is the 
analog of pairing#2 in Ref. [s^ replacing the basis func- 
tion (cos 

kx cos/cy), which represents nearest-neighbor 
pairing and transforms according to Aig , by cos k^ cos ky 
which corresponds to diagonal pairing and has the same 
symmetry. We will call this pairing if2\ It transforms 
according to Aig and it corresponds to intraorbital pair- 
ing. Following the previous steps we find that for pairing 
#2': 



MF 







^xy 


Afe 


\ 






^yy 





Afe 




Afe 





Cxx 


~£,xy 


V 





Afe 




~^yy I 



(B9) 



where Afe = V cosk^ cos ky, which is exactly the form of 
the s± pairing interaction proposed in Ref. [65l . 

In the base in which the tight binding Hamiltonian is 
diagonal we obtain: 
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where 



/ ei Ak\ 

62 Afc 

Afe -62 

V Afe -ei / 



(BIO) 



Note that this pairing operator corresponds to the s± 
pairing since it satisfies T^i(k) = — Vi(k + q) and there 
is no interband pairing. Thus, we have found that the 
s± pairing operator is possible in the two orbital model 
and transforms according to Aig. However, the Lanczos 
numerical calculations presented in the text suggest that 
in the region of physical interest the undoped ground 
state has to be connected to the ground state with two 
extra electrons via a pairing operator that transforms 
according to and, for this reason, the s± state is 
not favored. It only can prevail in the strong coupling 
regime of U and at the unphysical singular point U = U' , 
J = where the ground states in the doped and undoped 
regimes have the same symmetry and are connected by 
a pairing operator with Aig symmetry. 

For completeness, let us also consider a s± pairing op- 
erator frequently used in the literature.— i^i^i^ In this 
context it is assumed that A^ — — A^ = Aq which is 
independent of the momentum. Let us find whether this 
result is consistent with the symmetry of the two-orbital 
model. We start with ff^F ^^'^ working backwards the 
form of Hmf is found. By this procedure we obtain 



/ ei Ao \ 

£2 -Ao 

-Ao -£2 

V Ao -£i / 



(BU) 



that in terms of the original two orbitals corresponds to: 







A, 


Axy \ 






Axy 


Ay 


A, 


A^y 


^XX 


~Cxy 


Axy 


Ay 




~^yy / 



(B12) 



A^ = (u^ - t;^)Ao, 



2uvAo. 



(B13) 



(B14) 



(B15) 



We have found that uv transforms according to Big 
and — according to B2g; then this case corresponds 
to a linear combination of on-site intraorbital (#5) and 
interorbital (#8) pairings that transform according to 
two different irreducible representations of Z34h, i.e. Big 
and B2g. The coexistance of pairs with different sym- 
metries can occur only if the ground state with A'' parti- 
cles and/or the ground state with iV -I- 2 electrons are/is 
degenerate or nearly degenerate. Numerically, we have 
found that the ground states appear to be singlets and 
connected by an operator with symmetry B2g in the re- 
gion of physical relevance or Aig in the strong coupling 
limit. Only in the unphysical singular point J = and 
U = U' states transforming according to ^ig, Big, and 
B2g are very close to each other in energy. 

Summarizing, here it was shown that the s± pairing 
operator that pairs electrons along the diagonals of the Fe 
square lattice using the dxz and dyz orbitals does trans- 
form according to the Aig irreducible representation of 
the Z34h group that characterizes the symmetry of the Fe- 
As planes in the pnictides. Numerically we found that 
this is the pairing symmetry that prevails in the strong 
coupling region but that in the physical regime the pair- 
ing operator must transform according to B2g. In addi- 
tion, we have observed that the often-used momentum- 
independent approximation for the s± operator does not 
respect the symmetry of the Fe As planes. 
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